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Abstract. Linear processes are defined as a discrete-time convolution between a kernel 
and an infinite sequence of i.i.d. random variables. We modify this convolution by intro- 
ducing decimation, that is, by stretching time accordingly. We then establish central limit 
theorems for arrays of squares of such decimated processes. These theorems are used to 
obtain the asymptotic behavior of estimators of the spectral density at specific frequen- 
cies. Another application, treated elsewhere, concerns the estimation of the long-memory 
parameter in time-series, using wavelets. 
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f. roueff and m.s. taqqu 
1. Introduction 



Consider a linear process, that is, a weakly stationary sequence 

'^v{k-t)^k, keZ, where ^i^^(t) < oo 

and {^t, t £ Z} is a centered white noise sequence, that is an uncorrelated sequence with 
mean zero. We shall sometimes make the following additional assumptions on {^t, t G Z}. 

Assumptions A 

(A-1) {^t, t e Z} is a sequence of independent and identically distributed real-valued 

random variables such that E[^o] = 0, E[^q] = 1. 
(A-2) {^t, t € Z} satisfies and K4 =^ E[^^] - 3 is finite. 

We will allow decimation and consider, moreover, not one but N linear sequences, all using 
the same {^t, t G Z}. 

Definition 1. An array of N -dimensional decimated linear processes is a process admitting 
the following linear representation, 

^i,j,k = X] ^^'ji^jf^ -t)Ct , i = l,...,N, k ez, j >0 , (1) 

where {^j, t £ Z} is a centered weak white noise, (7j)j>o a diverging sequence of positive 
integers and, for all i = 1,... ,N and j > 0, {vij{t), t € Z} is real-valued and satisfies 

Remark 1. Zij^i^. involves three indices. The index i = 1, . . . ,N is used to define an A^- 
variate version, the index j labels the decimation factor 7^, and the index k corresponds to 
time. Because of the presence of the factor 7^ in ([T|), .^i.j,fc is not a usual convolution. It can 
be viewed as a decimated convolution of a white noise in the sense that, after convolution, 
one keeps only values spaced by 7j. A typical choice of decimation is 7j = 2-^, j > 0. 

Our goal is to study the asymptotic behavior of the sample mean square of Zij j., namely 
to find conditions on the kernels Vij , the decimation factor jj and normalization Uj , so that 
the normalized vector 



-1/2 



k=0 



converges to a multivariate normal Af{0, T) distribution. We want also to characterize the 
limiting covariance matrix T. Thus we are interested in the sum of squares of the Zij^k- 
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Such results are useful in estimation. In section [H for example, we apply our result to 
obtain a central limit theorem for the estimator fn(0) of the spectral density at the origin 
/(O) of a linear process. This CLT is compared to [4], Eq (3.9), as discussed in Remark [THl 
Another, more involved application, which involves wavelets, can be found in Q]. 

The paper is structured as follows. In Section [21 we indicate the main assumptions. The 
central limit theorems (Theorems [1] and [2]) for decimated sequences are stated in Section [3j 
Section [4] contains an application to the estimation of the value of spectral density at the 
origin (Theorem [3|) . Section [5] contains technical lemmas. Theorems [H [2] and [3] are proved 
in Sections EJ [7] and [8] respectively. 

2. Main assumptions 

Our assumptions will be expressed in terms of the Fourier series of the £^ sequences 
{vij{t), t € Z}, namely 

<,.(A) = (27r)-V2^^^_^.(i)e-iA*. (2) 

We suppose that for any i = 1, . . . , A^, as j — > 00, the Fourier series v* j concentrates around 
some frequency Aj^oo S [0, vr). By "concentrate", we mean that when adequately normalized, 
translated and rescaled around these frequencies, the series v* ^ converges as j — 00 to 
some limit functions v* with a uniform polynomial control (see Eq ([7]) and ^ below). 
Because of the particular structure of the 7^— decimation in ([T|), however, in order to derive 
the asymptotic behavior for the processes, we need to introduce sequences of frequencies 
(Ajj)j>o that satisfy some special conditions and converge to Xi^^o for alH = 1, . . . , A^. We 
shall first specify the conditions on the Fourier series v* j, the frequencies Aj^oo ^-nd the limit 
functions v^^, and then comment on these conditions. 

Condition C 

• There exist a A-dimensional array of frequencies {\ij)i£i^...,N,j>o valued in [0, vr) such 



that, for all i = 1, . . . , A, 

7jAjj G 27rZ+, for j large enough , (3) 

Aij Xi,oo , as j ^ 00 , (4) 

if Aj^oo = , then Ajj = 0, for j large enough . (5) 
and, for all 1 < -i < i' < A, 

if Aj^oo = Aj',00 , then Ajj- = Aj/j, for j large enough . (6) 
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• Moreover there exist J > 1/2 and a sequence of [— vr, 7r)-valued functions *&j(A) defined 
on A G M such that 

sup sup 7"'^'l<i(A)|(l + 7,|A - Xi,j\)' < oo , (7) 

i>0Ae[0,7r) 

nm77i/\*.(77iA + A,,,)e'*^W=<^(A) for ah AeR, (8) 
The fohowing remarks provide some insight into these conditions. 

Remark 2. Equations (jH) and dH) imply that the spectral density A i— > |'U*jP(A) of the 
undecimated stationary process 

concentrates, as j ^ oo around the frequency Aj^oo- In practical applications of the theorem, 
the limiting frequencies {Aj^oo, i € 1, . . . , N} are given. However, one can often easily find 
sequences (7j)j>o and (Ajj)j>o that satstisfy Conditions ^ and (jH) hold. In the particular 
case where the Aj^oo are such that q Aj^oo G 27rZ for alH = 1, . . . , and some positive integer 
q, one may take Xij = Aj.oo and 7j as a multiple of q. This happens for instance when the 
limiting frequencies are all at the origin, that is, Ai^oo = ••• = Aat^oo = and 7^ = 2^ . 

Remark 3. The presence of the phase function in ([5]) offers flexibility and implies that 
7j- '^l jilj ^ + -^0) converges to vl^{\) up to a common change of phase. Observe, 
however, that should not depend on i and thus, for > 1, Condition ([8]) is not equivalent 
to requiring that ^- jilj A + Ajj)| converges to |w*oo(A)| for all i. The presence of the 
phase $ is consistent with the fact that the asymptotic covariance matrix T defined in (|24p 
is invariant through a common phase translation of the functions v*^ for alH = 1, . . . , A^. 

Remark 4. Condition ^ states that if two limits Aj^oo and Aj/^00 are equal, then the Ajj- 
and Aj'j- which converge to them must coincinde for large enough j. Condition ([5]) has a 
similar interpretation. 

Remark 5. Conditions ([6]) and ([8]) imply that, for all 1 < i < i' < such that Aj^oo = Xi',oo^ 
lim77i[<.^](^7iA + A,,,) = Koo^](A) for ah AeM. (9) 
Here z denotes the conjugate of the complex z. 

Remark 6. Since Vij{t) is real valued, since Vij{t) is real- valued, we have 



<,(-A) = <,(A) 



(10) 
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Thus, Conditions ([7|) and Q imply that 

sup sup 77'^V*,.(A)|(l+7j||A|-A„-|)^<oo, 

j>0 Ae(-7r,7r) 



hm77l<,<,J(77'A-A,,,) = [^;*^<,,J(-A) for all AG 



j-»oo 



In particular, if Aj^oo = -^i',oo = 0, by ([5]), ([9|) and (fT2|) . we have 



[<oo<',J(A) = Koo<',J(-A). 



(11) 

(12) 
(13) 



Remark 7. Since (jj) is a diverging sequence and Xij — > Aj^oo ^ [0, vr), for any A G M, for j 
large enough, we have 'jJ'^X + Xij G [0,7r). Hence Conditions ([7]) and ([8]) imply that, for all 
i = l,...,N, 



sup|<^(A)| (1 + |A|)^< 



oo . 



(14) 



To better understand these assumptions, we start with a result on the asymptotic behav- 
ior of the cross-covariance function for the array ([T]). In this proposition, we set, without 
loss of generality, N = 2. 

Proposition 1. Let {Zi j^k, z = 1, 2, j > 0, G Z} he an array of 2-dimensional decimated 
linear processes as defined by Assume that Condition \(C)\ holds for some Xi^oo G [0,vr) 
and functions v* , i = 1,2, from M ^ Z. Then, for all k, k' G Z, as j — > oo, 



Cov (Zi,,- fc, Z2,j,k') ^ Ci,2 / wl,{X) e'^C^'-^) dA , 



where, for any i,i' G {1,2}, 



1 r- 



</(A) = ^ <oo(-A)<',oo(-A) + <oo(A)<,,oo(A) 



A G 



and 



if Ai,oo Ai' ,oo 

1 if Aj^oo = Aj'^oo = • 

2 if Aj^oo = Aj'^oo > 
Proof. Using ([1]) and Parseval's theorem, we have 

Cov {Zij^k, Z2j,k') = X] "ijiljk - t)v2j{-jjk' - t) 



(15) 



(16) 



(17) 



Kj^'2,i](A)e'^^^(*^'"^') dA 



(18) 
(19) 



We now consider separately the three cases Ai^oo 7^ A2,oo> Ai^oo = A2,oo > and Ai^ 
A2,oo = 0. 
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1) Suppose Ai^oo 7^ -^2,00- Then by ([7]), there is a constant C > such that 

rn 

|Cov [Zij^k,Z2j,k') \ < C-fj / (1 + 7j|A - Aij|)"''(l + 7j|A - A2, 

Jo 

^0 as i ^ 00 , (20) 



=,.,r'dA 



since jj — > 00, 5 > 1/2 and |Aij — A2j| has a positive hmit. 

2) Suppose Ai^oo = A2,oo > 0. Setting A = Jj'^^, + and using ([3]), we have 

rKK,]W e'^^-^^^'-'^) dA = ''"^ li'Kj^M'^ + Ai,,) e'^^'^'-'^) d^ 

Jo J~1j^i,, 

/oo 
Koo^](Oe^«('='-'^)de asi-00, (21) 
-00 

where the hmit follows from Conditions ([7]), Q and dominated convergence. 

Similarly we have 

fO roo 

[<,^](A)e'^^'^(^^'-'=)dA- / [^vloo]i-()e'^'^''-'^d( as j - 00, 

J —00 

by using ([H]) instead of ([9]). The last display, ([l9|) and yield 

/oo 
<2(6e'«(''-'Me asj^oo. (22) 
-00 

3) Suppose finally Ai^oo = A2,oo = 0. Setting A = gives 

e'^^-^(^'-'=) dA = ij'ivlAKlj'O e'^^''-'^ d^ 



/oo 
-00 



~oo 

by using Conditions dS]), ([7]), ([9]) and dominated convergence. The last display, (|20p 
and (IH]) yield i^. 

□ 

3. Main results 

We let — > denote the convergence in law. Our first result provides the asymptotic 
behavior of the sample mean square of an array of a decimated linear sequence under a 
global assumption on the behavior of the spectral density (the bound dZ])). A local version 
of this assumption is considered in Theorem [2j 

Theorem 1. Let {Zi j i^, i = 1,2, j > 0,k £ Z} be an array of N -dimensional decimated 



linear processes as defined by ([T]). Assume (A-2) and that is even for j large enough. For 
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each i = 1, . . . , N , we let Aj^oo denote a frequency in [0,7r) and v*^ a continuous 
function such that Condition \(C)\ holds. Then, for any diverging sequence (nj), 



1/2 



A;=0 



y2 



►Ar(o,r), 



(23) 



where T is the covariance matrix defined by 



^i i' = 47r Ciji 



dA, l<i,i'<N, 



where Cj^j' and w*-, are defined in (11) and \16\j. 



(24) 



Remark 8. From ()14p . it follows that the doubly infinite sum in (|24|) is well defined and 
bounded on A G M and hence F is well defined. 

Remark 9. The number of time indices k appearing in the centered sum in ()23p is Uj and 
asymptotic normalization occurs as j and Uj tend to oo. 

Remark 10. The presence of the factor 7j in ([T]), and hence of decimation, is essential for the 
Central Limit Theorem to hold in this generality because it ensures that the dependence 
of the Zj j^fc's decreases sufficiently fast as j oo. Decimation of this type is typically 
encountered in settings involving wavelets, or more generally filter banks, see Q]. 

Remark 11. In applications, the expectations in (123]) . which depend on j, will be approx- 
imated by quantities that are independent of j. To see why this is possible, observe that. 



applying Relation (jlSp in Proposition [T] with k = k' = and i 



1, . . . , N, we get 



hm E o] 



|<^(A)|'dA<oo. 



(25) 



Thus, when the convergence rate to this limit is fast enough, the expectations in (j23p can 

2 



be replaced by <oo(^) 



dA, i = I, . . . , N, which does not depend on j. 



We have assumed in ([7]) a bound for v* j{X) for A € (— vr, vr). This bound implies that the 
spectral density of the process Zij^, defined in ([T|) is bounded on (— 7r,7r). We shall weaken 
this assumption by only assuming a local bound around the frequency Xij as follows. 

Theorem 2. Assume that all the conditions of TheoremU\ hold except that ^ is replaced 
by 

sup sup 7~'/V*^.(A)|(l + 7j|A-A,j|)'5 <oo, (26) 

i>0 |A-A,,^|<e 
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where e > is arbitrary small. Suppose in addition that 

n/^ ^ 1(|A - Ai,oo| > e) |<j(A)P dA ^ as j ^ oo . (27) 
Then the conclusion of TheoremUl that is, the CLT Il2cl\} . still holds. 



Remark 12. Since ([7]) is replaced by the local condition (|27p . we impose the additional 
condition (j27|) on the growth of Uj. This condition does not appear in the conditions of 
Theorem [H where it was only required that nj ^ oo. 

4. Application to spectral density estimation 
Let {X{u), M € Z} be a standard linear process, 

X, = ^a(u-t)et , (28) 

tez 

where {^t, i G Z} is a centered weak white noise with unit variance and {a{t), t G Z} is 
real-valued sequence such that Ylk ^ ^ with Fourier series 

a*(A) = (27r)-i/2 ^ ^-ixt _ (39) 

Then {Xk, S Z} admits the following spectral density 

/(A) = |a*(A)p, Ae(-7r,7r). 



For simplicity, as in Section 3 of [2], we consider the problem of estimating /(O) from 
observations Xi, . . . , Xn • 

Let us denote hy W a. bounded M ^ M function with compact support and by W its 
Fourier transform, 



/oo 
W{t)e-'^^ dt . 
-00 



Let (7j) be any diverging sequence of even integers. 

We let N = 1, and Xij = Ai,oo = for all j > 0, which yields (g]), §i) and ([6]) in 
Condition [(C)| 
Define 

We assume that 
(H-1) As A ^ 0, 

/(A) = /(0) + O(|Ap). (31) 
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(H-2) the support of W is included in [-1,0], sup^gjj + ICI)^ < oo with /? > 1 

and 



W{X) 



dX = 1 



(32) 



Assumptions (H-1) and (H-2) are related to the standard bias control of kernel estimates 
of the spectral density (see Lemma [8] below). 
Define 

n, = h-\n + l)] . (33) 
For all /c = 0, rij — 1, since W{k — "ij^u) vanishes for u < and n > n + 1, we have 



Z 



-1/2 



In other words, {Zij^, k = 0, ...,nj — 1} can be computed from the n observations 
Xi, . . . ,X„. Thus 



/n(0) =nTi J] 



fc=0 



can be used as an estimator of /(O). The following theorem provides a central limit result 
for /„(0). 



Theorem 3. Assume (H-1) and (H-2) with (3 > 2. Let (7^) be a diverging sequence of even 
integers such that '^n'^n — > 00 Then, as n — > 00, 



E 



/„(0)1 =E[z2„,o] =/(0) + O(7n')- 



If moreover (A-2) in Section\M holds and 



,1/2^1/2-2/. 



then 



where 



(34) 
(35) 
(36) 



a2 = 27r/(0)2 / iY,\wiX + 2p7T) 



dX 



Remark 13. Our CLT (j36p can be compared with [j, Eq. 3.9], although the estimators 
are different since ours involve a decimation and the one m a is express. a wei ffhted 
integral of the standard periodogram. In (I36p . our 7„ has a role similar to the q = Qn 
for their estimator. Our bias estimate ()34p has a faster decrease than the corresponding 
one 0(g~J_) in 0], see the last display in their Section 3. Our conditions also differ from 
those of [3]- Our conditions on the weight sequence a{t) is much more general, since we 
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assume a polynomial decrease neither of this sequence nor of a{t) — a{t + 1), as assumed 



for the corresponding (causal) sequence (^j) in Assumption 2.1 of [4]. Standard results on 
spectral estimation (see e.g. l|, Theorem 9.4.1]) usually assume the even stronger condition 
\a{t)\ < 00. On the other hand we do assume that the noise sequence {^t} has fourth 



finite moment (see (A-2) ) while only a finite 2 + /? moment (with (3 > arbitrary small) is 
assumed in Q]. It is an open question whether similar moment condition can be used for 
our estimator /n(0). 

5. Technical lemmas 
The following lemma will be used several times. 

Lemma 1. Let g be a {2tt) -periodic locally integrable function. Then for all positive integer 
7, the function defined by 

7-1 

g^(A) = J]5(7-^A + 2pvr)) 

p=0 

is (2tt) -periodic. Moreover, one has 

r giX)dX = r' r 57(A) dA. (37) 

J —n J —TX 

Proof. Observe that, for all A G M, 

7-1 

g^{X + 2^) = Y,9 (7^' (A + (p + 1) 27r)) 

p=0 

7-1 

= Y.9 (7"' (A + 2p7r)) + g (7-^ (A + 72^)) = g^[X) , 
p=i 

since g{'y~^ {X + 'y27r)) = g{'y~^X) by the (27r)-periodicity of g. Hence g-y{X) is (27r)-periodic. 
With a change of variable, one gets 

r-27r /.27r7-l /■27r 



7 



1 / 57(A) dA = V / 5(e + 2P7TJ-') de = / g{0 d^ . 
Jo Jo Jo 



Relation ([37|) follows by (27r)-periodicity of the integrands. □ 

The next lemma relates the rates of decrease of two functions with the rate of decrease 
of their convolution. 

Lemma 2. Let 5 > 0. For all T > and to G [0,T/2), we let hT,to{t), t e M be the even 
and T-periodic function such that 

hT,toit) = {l + \t- to\)-^ for all t G [0,T/2] . 
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Let g be an integrable non-negative function on M such that 

9{t) < co\t\~^~^ for \t\ > 1. (38) 

Then there exists a positive constant c, depending only on 5, \\g\\i = f^^g{s) ds and cq 
such that, for all T > 0, to £ [0,T/2), and t G [0,T/2], 

/oo 
g{t - u)hT,to{u) du < c{l + \t - to\)~^ . (39) 
-oo 

Proof Let to e [0,T/2). We shall use the bound, vahd for all t G M, 

9 * hT,to{t) < \\g\\i . (40) 

This bound yields (I39p only for t close enough to to • We shall derive a different bound valid 
only for t G [0, T/2] with \t — to| ^ 2, namely 

g * hT,to{t) < 2"-^ Co 5-^ \t - to\-^ + llslli (1 + |t - tol/2)-'' • (41) 

Applying (gOD for |t - to| < 2 and gH) for |t - to| > 2 yields (l39l) . 

Hence it only remains to establish ([3TD for t G [0,r/2] with \t - to| > 2. We shall 
suppose that t G [to + 2, T/2] (the case t G [0, to — 2] is obtained similarly). Let u such that 
\t-u\<\t- tol/2. Then t - (t - to)/2 < u < t + (t - to)/2 and thus, using that t < T/2 
implies t + t/2 < T - t/2, we get (t + to)/2 < n < T - (t + to)/2. Observe that the middle 
point between (t + to)/2 and T — (t + to)/2 is T/2. Since hT,to{u) is decreasing on [to, T/2], 
and symmetric around T/2, we get hT,to{u) < hTfy^{{t + to)/2) = (1 + |t — to|/2)~'^. Hence 
we may bound hT,to{u) by 1 for |t — m| > |t — to|/2 and by (1 + |t — to|/2)~'' otherwise, which 
gives 

/oo r 
g{t - u)hT,toiu) du < / g(s)ds + ||(7||i (l + |t-to|/2)-'5 . 

-oo J\s\>\t-to\/2 

Since |t — to| > 2 we may apply the bound ([38]) in the integral of the RHS of the previous 
display. Hence we get (j4ip . which concludes the proof. □ 

The following lemma is used, in particular, to bound fj in the proof of Theorem [2j It 
will be used again in the proof of Lemma [6] below. Applying it, one can bound gj^^.{\) 
independently of j and A, where gj^^. is defined as in Lemma [J with g replaced by gj, and 
the sequence gj satisfies a uniform bound of the form ([7]), namely 

sup sup |gj(A)|(l+7||A| -Ajl)*^ < oo , 

i>0 Ae[-7r,7r) 

with \j Aoo G [0, vr) as j — > oo. 
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Lemma 3. Let 6 > 1/2. Then 



sup sup ^(1 + ||t + 2p^| -i'D^^-^ < oo . (42) 



Moreover, as it —> oo, 

sup sup V(l + \\t + 2pTT\ - t'\)-\l + \\t + 2pTT\ - t' -u\)-^ ^0 . (43) 

Proof. Let S{t,t') = Epez(3^ + II* + 2pvr| - t'\)~'^^. Since, for any t', t ^ S{t,t') is (27r)- 
periodic we have 

sup sup S{t, t') = sup sup S{t, t') . (44) 

tGMi'GM t'eM t:\t~t' \<-K 

Suppose that t, e M are such that |t — 1'| < vr. Then for any a, 6 E M, we have, if t + a > 0, 

||t + a| > |a-6| -vr , (45) 

and, if t + a < 0, 

\\t + a\-t' -h\ = \t-t' + 2t' + a + h\> \2t' + a + b\ - n . (46) 

Adding Svr to each of the last two displays with a = 2pTr and 6 = 0, we get that, for all 
\t -t'\ <TT and p G Z, 

(37r+ Ijt + 2p^| -t'\)-^^ < (2|p|7r + 27r)-2'5 + {\2t' + 2p7r\ + 2tt)-^^ . (47) 



Since Epez(2|pK + ^)~^^ < °^ and supt,giR ^pg^(|2i' + 2p-K\ + 7r)~2<5 < qo. Relation (|44 
gives that sup^g^ sup^/gjg S{t, t') < oo and follows. 
We now prove (03]). Let 



S{t,t',u) = ^(37r + ||t + 2p7r| - t'|)~'^(37r + \ \t + 2pTr\ - t' 



As above, we have 

sup sup 5j(t, t') = sup sup Sj{t,t') . (48) 

teKt'GR t'eM t:\t-t' \<TT 

Suppose that t, t' G M are such that \t — t'\ < vr. Adding Sir to (05]) and ([I6]) with a = 2p7r 
and b = u, we have 

(37r + ||t + 2p7r| - t' - li])"^ < (|2p7r - u\ + 2tt)~^ + (|2t' + 2p7r + n| + 27r)~^ . 
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Using ([18|) . (jT7|) and the previous display, we obtain 

sup sup S{t, t\ u) < ^(|2p7r| + 27r)-^(|2p7r - u\ + lir)-^ 



pGZ 

) 7 y\^b Z,p/i I -|- Z7I ^ -|- ^p/i -\- U\ -f- ^71 ) ^ 



+ sup ^(|2t' + 2p7r| + 27r)"'^(|2i' + 2p7r + n| + 27r)" 
+ sup ^(|2p7r| + 27r)~'^(j2t' + 2p7r + li] + 27r)-'^ 
+ sup^(|2t' + 2p7rj +27r)~^(j2p7r-'u| + 27r)~'^ . 



Since the three functions in t' appearing in the right-hand side of the last display are vr- 
periodic the supj/gjj can be replaced by sup|^/|<^/2- Since \2t' + 2p7r| > \2p7r\ — vr and 
\2t' + 2p7r + n| > \2pTr + u\ — it for \t'\ < 7r/2, we thus obtain 

sup sup as u —>■ oo , 

tm i'GM ^ 

which conclude the proof. □ 

The following lemma will be used in the proof of Lemma [H 

Lemma 4. Letp be a positive integer. For all C'P -valued function g € L^(— vr, vr) and n > 1, 
define 

where = {2tt) J^^ ^^^'^ dA and \ ■ \ denotes the Euclidean norm in any dimension. 
Then, for all gi and g2 in vr, vr), 

1 /2 

|M„(gi) - M„(g2)| < (^y" |gi(A) - g2(A)|2 dA^ . (50) 

Moreover, for all g in L'^{—-k, vr), as n —>■ oo, 

/ rn \i/2 
M„(g)^ / |g(A)|2dA . (51) 



Proof. See p. Lemma 1 (Appendix B)]. □ 

The following lemmas are used to compute the limiting covariances ()106p and (j24p . 

Lemma 5. Let {Zi j ^, i = 1,2, j > 0,k £ Z} be an array of 2 -dimensional decimated 
linear processes as defined by ([T]). Assume (A-2), Then for all j > and all n > 1, one has 

-Gov ( ^ Zl^^„ J2 Zlj,k I = 2A,(n) + K^B,{n) , (52) 



vfc=0 k=0 
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n-1 



r=— n+1 



n-1 



weZ T=-n+l 

Proof. Using a standard formula for cumulants of products, we have 



(53) 



(54) 



n-1 



n—1 n—1 



\k=0 



k=0 



J^^Cov(Z2j.fc,Z|j.fc,) 
fc=Ofc'=0 
n—1 n—1 



fc=0 k'=0 
n—1 n—1 

+ ^ ^ cum (^lj,fc, -Z'lj^fc, -^2J,fc', ^2,j,k') ■ 
k=Ok'=0 

By definition of {Zij^i^, i = 1, 2, € Z} the covariance and the fourth-order cumulant in the 
previous display read respectively 

Gov {Zij^k, Z2,j,k') = 'Y'l^ijiljk - t)v2ji^jk' - t) 



and 



cum {Zij^k, Zij^k, Z2,j,k',Z2j,k') = ^4 ^ vl j{-fjk - t)vl j{'^jk' - t) . 
The two last displays thus give ([52]) . 



□ 

Lemma 6. Let Aj{n) and Bj{n) be defined by i53\) and |5-^p , respectively, and v* ^ by f^. 
Then the following inequalities hold for all j > and all n > 1 : 

2 



Aj{n) < 2-K 



7,-1 

77'EK/<,](77nA + 2vrp)) 

7i-l , 



dA 



p=0 



dA . 



(55) 



(56) 



Suppose moreover that 'jj — > oo as j ^ oo, -jj is an even integer for j large enough and that 
Condition \(C)\ holds for some Aj^oo £ [0)^) ^.''^d, M — > Z functions v*^, i = 1,2. Then, as 
{n,j) (oo,cx)), 

2 



Aj{n) ^ 27r Ci,2 [ 



J^^^2(A + 27rp) 



dA , 



(57) 
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where Ci^2 '^^'d 0,'^^ defined in (11) and M0\j respectively. Moreover, 



lim sup|Sj(n)| = . (58) 

Proof. Step 1. Using properties of the convolution of square summable sequences, we have, 
for all t G Z, 

../-^ J —IT 



For any r G Z, applying Lemma [T] to the (27r)-periodic function A ^ j{X)v2 e^'^^^'^, 
one gets 



E 

where 



i,A^)=l-\X + 2TTp). (59) 
Using the notation of Lemma HI we can express Aj (n) defined in (|53|) , as 

A,{n) = {Mn{g,)f (60) 

where 

7i-l 

5,(A) = (2^)V2^-i <,(^,,p(A))^;2*,,.(0,p(A)) . (61) 

p=0 



The bound (I50p in Lemma |4] with gi = (7^ and g2 = thus gives ([55 
Step 2. Let us now show (I56p . Since 



^2,,(7.-r + n) = (27r)-i/2 / ^;;*,,(A) e^^^^^'^+^^dA 



we can apply Lemma [T] to the (27r)-periodic function A 1-^ 'I'l j(A) e' '^'^'''j for all u and r 
in Z and get 



V2Aljr + u) = {2nr'/' [ 7^' E ^2,,(0.p(A))e'<-W j e^^dA 
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Using the Parseval formula, we get, for all ti G Z, 



n-l 



^ (1 - \T\/n) 4,-(7jT + u)<Y, + 



u] 



T=— n+l 



7,-1 
p=0 

^ j\{^j''Y.H,MU^))\^ dA. 

Observing that the resulting bound is independent of u S Z and using the Parseval formula 
X^ugz^i j(^) ~ /-TT '^1 jW dA, we obtain the bound ([56]) for Bj{n) defined in ([M]) . 
Step 3. We now establish the limit ([57]) successively in the cases Ai^oo 7^ ^2,00 and Ai^oo = 
^2,00- The (27r)-periodicity oivljV2j and Lemma [T] entail that gj (defined in ()6ip ) is (27r)- 
periodic. By definition of M„ in Lemma HI it follows that, for any j > and any r € M, 



M„((7j) =M„(5f)) with g'f\x)=gj{X-T), A 



(62) 



since the modulus of the Fourier coefficients of gj and g^ are equal. In the following we 
will take r = tt^j. Observe that, for all p G {0, . . . , 7^ — 1}, A € (0, 27r) and j > 0, 

^j- p(A - 7r7j) G (-TT, vr) . (63) 

Consider the case where Ai^oo 7^ ^2,00, which, by (jH), implies 

7j|Aij — A2j| ^00 as j ^ 00 . (64) 

Using dMD, dSI]), dSl]), dM]) and (dH), we have, for some constant C > 0, for all j > 0, 



sup 

Ae(0,27r) 



^(-^^^(A) = sup \gj{X-7T^,)\ 

AG(0,27r) 
7j-l 1 

< C sup ^(1 + \\t + 2TTp\ -t'\)-^ 



X (1 + ||t + 27rp| - i' - 7j|Aij - A2j||)~* 
— as j — > 00, 

by ([M]) and (fi3|) in Lemma [3l Applying ([60]) . (f62|) and the bound ([50|) in Lemma [3] with 
gi = gj and g2 = yields 



A,{n) = {AUgf^^^)f < I \g) 



,{'^7j)|2 



dA . 
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The last two displays and the (27r)-periodicity of gj imply Aj{n) ^ as j — > oo. This 
proves ([57|) since by p!7|) . Ci^2 = when Ai^oo ^2,oo- 

We now consider the case Ai^oo = ^2,00- By Condition we have Xi j = A2J for all j 
large enough. Let pj = 7jAij/(27r) +jj/2 so that 

Ai,j = A2,j = 27rj~^pj - vr . (65) 

By Condition ^ and since 7j is even for j large enough by assumption, we get that pj is 
an integer for j large enough. Writing 

7i-l [7,/2]-l 7,-1 [7j/2]-{7,-p,)-l -yj-Pj-l 

1;= E + E = E + E . 

p=0 p=0 P=[7j72] g=-{7j-pj) r=[7j/2]-pj 

where q = p— {jj —Pj) and r = p—pj and observe that, with these definitions, ()59p and (165^ . 
we have = 7j"^(A + 27rg) = 7j"^[(A - 7r7j) + 27rp] + 27r7j"Vj - vr = Cj,p(A - 717^) + Aij, 

and, similarly, Cj,gW = ?j,p('^ ~ ''^Ij) ~ so that g^^"'^\x) = gj{X — ir^j) defined in ([62]) 



and (j6T]) . can be expressed as 



5f^^-)(A) = (2.)^^ 



[7j/2]-{7i-Pi)-l „,* — T" 
^lj^2,i 



E ^fe.(A)-Ai,,) 



'?=-{7i-Pi) 



r=[7,/2]-p, 



(66) 



Since limj_^oo 7j = 00 and, by Condition ([!]), limj^oo Aj j = Aj_oo G [0, vr), we have limj^oo 7j ^Pj G 
[1/2, 1) and thus 

— (7j — Pj) — > —00 and — pj — 1 ^ 00 , (67) 

namely, that in ()66p . the upper limit of the first sum tends to 00 and the bottom limit of 
the second sum tends to —00. We now consider the remaining limits. If Ai^oo = A2,oo > 0, 
then limj^oo 7j^^Pj falls in the open interval (1/2, 1) and thus 

[7i/2] - (7i - Pj) -I ^00 and [7^/2] - pj ^ -00, . (68) 

If Ai^oo = A2,oo = 0, using dS]), ([65]) implies pj = 7j/2 and thus, for j large enough so that 
Pj is integer-valued and 7j even, 

[7i/2]-(7i-Pi)-l = -l and [7^/2] - p,- = . (69) 



In view of ([SSD, Conditions ([6]), (UD and ^ (which imply (l9l). (fTT]) and (fT2]) ). (|671) . (|68]) . (|69]) 
and dominated convergence yield, for all A G (0,27r), 

5f^^\A)^goo(A) asj^oo, (70) 
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where 

5oo(A) = (27r)i/2 

if Ai,oo = A2,oo > 0, and 
5oo(A) = (27r)i/2 



EKoo4oo](-A - 27Tq) + 5]K,oo^^2,oo](A + 27rr) 



E Koo^2,oo](-A - 27rg) + EKoo^.ooKA + 27rr) 



.g=— oo 



r=0 



if Ai^oo = A2,oo = 0. By definition of w* in and using ([13]) . one has \v\^v\^{—X) 
\P\ oo^2 oo](A) = 2(A)) the two previous displays read 



500(A) = (2^)1/2 Ci,2 E<2(A + 2^p) 



(71) 



Conditions ^ and ([8|) imply (fH|) . thus that goo(A) is bounded, hence square integrable on 
A E (— 7r,7r). Moreover, applying the same dominated argument as above, one has 



lim 



dA = 



(72) 



(73) 



One gets bv (pD . (I50]) and (I72D 

/TT 
l5i(A)-5oo(A)|^ dA^O asj^oo. 
-TT 

By applying the limit (f5TI) with g = gioo, one gets 

Mn{g^f ^ r boo(A)pdA (74) 



as n ^ 00. Hence, setting Mn{gj) = {Mn{gj) — Mn{goo)) + Mn{goo), the limit (j57|) follows 
from dm]), jnl), dZH]) and jTlD. 

Step 4. We now establish the limit (jSSp . By Condition ([7]), we have 

sup / 'Wi,(A) dA < CX3 . (75) 

i>0 J~7T 

Using similar arguments as above and Condition d?]), we have 

2 



/ 77^ E K(e,.p(A))| dA = 7-W E 



7"'^'k2%(e,,p(A-7r7,))| I dA 



r-27r n^-l 



^^77'^"f E(l + ll^ + 2^(^'-7.-/2)|-7jA2,il)"M dA (76) 



Using that ||a + 6| — c| > — c| — \a\ and A € [0, 27r], we have 

||A + 27r(p - 7j/2)| - 7jA2j| > |27rjp - 7j/2| - 7jA2j| - 27r 
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Take j large enough so that 7j is even. Since A2J G [0, vr), as p € {0, ... ,7^ — 1}, 27r|p — 
7j/2| — 7jA2j is a sequence of numbers with lag 2tt and belonging to [— 7j7r,7j7r] and can 
thus be written as a sequence 27rg + c, where q belongs to {—7^/2, . . . , 7j/2} and c to [— vr, vr] 
so that 

\2tt\p — 7j/2| — 7jA2,j| — 27r > 27r|gj — Svr . 
From the last two displays, we have 

(5^ + ||A + 2n{p - 7,/2)| - 7,A2,,|)-' < (2^)-'(l + kl)"' , 

with q describing {—7^/2, . . . ,7j/2} as p describes {0, . . . ,7j — 1}. Inserting this bound 
in ([76|) . we get 

•^-'^ \ P=0 / \q=0 

for some constant C not depending on j > 0. Since the last right-hand side of the previous 
display tends to as j ^ 00 for any 5 > 1/2, with (f75]) and ([56]) . we obtain ([58]) . □ 



Remark 14. The factor 47r = 2 x 27r in (j24p is due to the factor 2 in the right-hand side 
of ([52|) and the presence of 27r in the right-hand side of ([57]) . 



Corollary 1. Let {Zij^^, i = 1, 2, j > 0, /c G Z} be an array of 2- dimensional decimated 
linear processes as defined by HP- Assume (A-2), that is even for j large enough and 
that Condition \(C)\ holds for some Xi^oo G [0, vr) and M ^ Z functions v*^, i = 1,2. Then, 
for all k, k' € Z, as j 00, 



Cov(z2,.„Z2V) -2C?,2 (y"^<2(A)e'^(^'-'=)dA)' , 
where Ci^2 and '^'^e defined in [T7\ ) and fT6\). and, as {n,j) — > (00,00), 



/n-l n-l \ 

-Gov E4.>'E4.> -4vrC,,2 / 
^ \k=o k=o I •'-^ 



E<2(A + 2vrp) 



dA 



(77) 



(78) 



Proof. Setting n = 1 in ([52]) and ([53j) and replacing vij{—t) by vij^'jjk — t) and V2j{—t) 
by V2j{'~fjk' — t) so that .Z'lj-^o is replaced by ■Z'lj-^fc and .^2,j,o by -^2,j,A:'; we get 

Gov (Z2^. „ = 2Aj{l) + K4i?,(l) 



2 ( E ^ij(7j^ - *)^^2,i(7jA:' -t)] + K45i(l) , 
Vtez / 



and thus ([77]) follows from (pD . (flSI) and (fT5[) . 

Relation (178]) is obtained by applying Lemmas [5] and El 



□ 
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6. Proof of Theorem [T] 

We first establish two Central Limit Theorems which will be used in the proof of Theo- 
rem [TJ The first involves a sequence of linear filters of the sequence {^t, t € Z}. 

Proposition 2. Define, for alii = I, . . . , N , j > 0, k G Z, 

Zi,j = Y.Vi,,{t)Ct, (79) 

tez 

where for all i = I, . . . , N and j > 0, {vij{t), t ^ Z} is real-valued and satisfies Y^t& "^fji^) < 
oo and {^t, t € Z} satisfies (A-1). Assume that 



lim sup = for all i = 1 



,...,N, 



lim = fof 0,11 i,i' = 1, . . . , N 



(80) 
(81) 



where is a N x N given matrix. Then, as j ^ oo. 



Z 



N,j 



^AA(0,S) 



(82) 



Remark 15. A study on the_weak convergence of such sequence without assuming Assump- 
tion 



(A-1) can be found in 



Proof. This is a standard application of the Lindeb erg- Feller Central Limit Theorem. Using 
the Cramer- Wold device for vectorial Central Limit Theorem and since for any column 
vector C G M^, the linear combination "^j^i Ci^i.j can be written as in (I79p with Vij(t) 
replaced by Vj{t; C) = J2iLi d^ijit), which satisfies 

lim sup \vj{t; C)\ = for all i = 1, . . . , N , 
liiny2iv,{t;C)vj{t;C)) = C^^C , 



it is sufficient by ()80p and (jSip to prove the result for = 1, in which case we simply denote 
vij{t) by Vj{t) and S by a^. Let {nij) be a sequence of integers tending to infinity with j, 
such that 

\t\>mj 

We now show that the Lindeberg conditions hold for the sequence X]|t|<mj '^ji't)^t- The first 
holds because, by ([HT]) . 

lim y Var(?;, (t)et) = lim V ^;2(t) . 



t|<mj 



j-»oo 



t|<mj 
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The second holds because, for all e > 0, 

E [v]m^ i{\v,mt\ > ^)] < ( E^K*) ) ^ fe' i(ieoi > e/s,)] , 

\t\<rnj VteZ / 

where Sj = supjg^ bj(OI) ^^^j by ([80|) . the right-hand side of the last display tends to as 
j ^ oo. This concludes the proof. □ 

The second Central Limit Theorem deals with m-dependent arrays. Recall that {Yk} is 
said to be m-dependent if, for all p > 1 and all ki, . . . ,kp such that ki + m < k2, ■ ■ ■ , kp-i + 
m < kp, Ykj^, . . . , Yfcp are independent. 

Proposition 3. Let m be a fixed integer and (rij) a sequence of integers such that nj — > oo 
as j oo. Let {Yj^^, k = 0, . . . ,nj — 1, j > 1} be an array of W^-valued random vectors, 
such that, for each j > 0, {5^',fc, fc = 0, . . . , nj — 1} has zero mean, is strictly stationary and 
m-dependent. Assume that there exists a dxd matrix T and a centered stationary -valued 
process {Y^, > 0} with finite variance such that 

Yj,. ^ Y. as j ^oo , (83) 
lim Gov (Yj,k, Yj^k') = Gov (Y^, Yfe/) for all k,k' >0 (84) 

lim lim Gov f r = T . (85) 

Then we have, as j ^ oo, 

nj — 1 

n7'/'5^y,-,fc^AA(0,r). (86) 

k=0 

Proof. We may suppose that d = 1 since the vector case follows by the Gramer-Wold device. 
For convenience, we set F = o"^. Let s be a positive integer larger than m. We decompose 
of YlkLo^ ^j,k ill sums of random variables spaced by m, as follows: 

rij-l pj pj 

k=0 k=0 k=0 

where 

s-l m-l Qj-l 

i=0 j=0 i=0 

and where pj and qj are the non-negative integers defined by the Euclidean division 
Uj = pj{m, -|- s) -|- qj with gj € {0, . . . , m -|- s — 1} . 
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The m-dependence and the strict stationarity of the sequences Yj^, ensure that for all 
J > and all s > m, S'j^_^ and tJ^'' are sequences of centered independent and identically 



distributed random variables. Hence, by 



we have 



J— >CX3 



k=0 



lim 5; Var (pj'^'s^l) = liiTi Var (5];^) = Var (s(^) 



where 



i=0 



lim Var ( pj Tj'^^ 



k=0 



lim Var ( T, 



Var ( S^'") 



and 



lim sup Var (^R^^^ < max |Var (^S^*)^ : t = 0, 1, . . . , m + s - l| < 00 . 

In addition, for any e > 0, 

Pj 



fc=0 

and hence, since pj ^ 00, 

Pi 

hmsup^E (pTV25(J)2j(^-i/2|^(s)| ^ ^) 



E 



{S^iflipJ"'\S^l\>e 



l/2|r.(s)| 



fc=0 



But using (j83|) . we have 



< inf lim sup E 

A/>0 J ^00 



5]J^SW as j^oo 



(5]5)21(|S]5|>M) 



(89) 



(90) 
(91) 



(92) 



Hence, denoting by (pM some continuous M+ [0, 1] function satisfying l(x < M/2) < 
(pMi^) < t{x < M), so that x^4>m{x) is continuous and bounded, we have 



E 



> E 
^ E 
^ E 



(4:J)Vm(5(:oV 

(S(^))2<^m(sW) 



as J ^ 00 



(SW)2 



as M — > 00 . 



Using (IMI), we have E \{S^'of 



E [(S'-'^-')^] as j ^ 00 and hence, for any M > 0, 



lim sup E 



(5«)^1(|S]5|>M) 



E 



(S(«))2 



lim inf E 

j-»oo 



(S]5)'1(|5]J|<M) 
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The two last displays and (j92p imply the second Lindeberg Condition, namely, 
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limsupJ^E \{pf^sflfl{pf'\sf^\ > 



for any e > . 



Using this and (j89p . we may apply the Lindeberg- Feller CLT for arrays of independent r.v.'s 
and we obtain 



1/2 ■ 



5:5(5 ^AA(0,Var(s 



k=0 



i(s) 



lim s^^Var ( S^'^ 



By (l85|l . we have 



Using ()93p and ()88p with j ^ 00 and then ()94p with s — > 00 yields 

Pj 



(93) 



(94) 



n. 



-1/2 



ES]S (»■(■" + »)-'Var (S«-))) _^^M(<,y) 



k=0 



On the other hand, by ([87]), ([90]) and §T^, we have 



lim sup E 



''3 



1/2 



< (?Ti + 5)-^Var (^S^™)) ^ 



fc=0 k=0 

as s ^ 00. Using the last two displays and Theorem 3.2], we obtain (j86p . which concludes 
the proof. □ 

Proof of TheoremUl The proof is in three steps. We show in a first step the convergence of 
the process [^ij-,fc, i = 1, ■ ■ ■ , N, k = 0, ... ,nj] as j ^ 00 towards a Gaussian limit. In the 
second step we prove Theorem [1] under the additional assumption that the sequence 



def 



, k = 0, . . . ,nj — 1 



is m-dependent. The third step exhibits an m-dependent approximation and extends the 
m-dependent case to the general case. The proof uses a number of auxiliary results proved 
in Section [5l 

Step 1. We shall apply Proposition [2l By Relation (jlSp in Proposition [H we get, for all 
= 1, . . . ,N and all k, k' € Z, as j — > 00, 



Gov {Zij^k, Zi'j^k') Cj^j' 



w*i,{X) e'^('='-'=) dA 



(95) 
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Moreover, by d?]), one has, for all i = 1, . . . , A^, 

svip\vij{t)\ < 2 (27r)-i/2 r jx)\ dX 
tez Jo 

7j{7r-Aij) 



< 2 (2vr)-i/S7V2 



(1 + |A|)-MA, 



which, by tends to as j — > oo for any 6 > 1/2. Hence, by Proposition [2l for any p > 1, 
any ii, . . . ,ip G {I, . . . , N} and any ki, . . . , kp G Z, we have, as j oo. 



►AA(0,S) 



where S is the covariance matrix with entries given for all 1 < n,n' < p hy 



W, 



' , ,(A) e'^('="'-'=")dA . 



Expressing this integral as Y^pH^+^T^^ S^* 



(96) 



j;<,^,(A + 2pvr) e^^(^"'-'=")dA. 



The convergence (j96p can be written equivalently as 



N,j,. 



z. , 



(97) 



as j — > oo, where {Z^ = [Zi^oo,k ■ ■ ■ Z]\f^ao,k]'^ , A; > 0} is a stationary Gaussian M^-valued 
process with spectral density matrix function D with entries 

A/ (A) = J2<>^J^ + 2^^^)' 1 < < ^ • 

Step 2. In this step, we prove ()23p . assuming that for each j > 1, {Yj k, k = 0, . . . ,nj — 1} 
is m-dependent. We shall apply Proposition [3] under this additional assumption. We thus 
need to show that ([83]) - ([85]) hold. Relations (p5|) . ([97|) and the continuous mapping theorem 
imply (j83p with 



l,oo,fc 



E[Zl 



oo,fcJ 
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Since Z. is Gaussian, we have, for all k,k' > and all = 1, . . . , N , 

.2 



Gov (Yfc^i, Yk'^i') = 2Cov {Zi^oo,k, ^i',oo,fc') 



2 



Hence Relation (j77p in Corollary [T] and the previous display yield (|84p . The final con- 
dition (|85p follows from Relation (j78p of Corollary [1] with a covariance matrix F with 
entries (|24p . Applying Proposition [3] then yields (j23p . with F given by 



Step 3. Let K{t) be a non-negative infinitely differentiable function defined on t G M whose 
support is included in [—1/2, 1/2] and such that A'(0) = 1. We will denote by K its Fourier 
transform, 

/oo 
K{t)e-'^^ dt . 
-oo 

Observe that, by the assumptions on K, -fC(C) decreases faster than any polynomial as 
1^1 — > oo. In particular i^(^) is integrable on ^ G M and, for all i G M, 

1 f°° ^ 

m = ^j i^(Oe'«*de. (98) 

The function K will be used to approximate the v^"^^ sequence by a sequence whose 
dependence structure can be controlled. We thus define, for any i = 1, . . . , N, j > and 
m > 1, 

which vanishes for all \t\ > and we define Z^^^^ and v* accordingly. 

Let m be a fixed integer. Then, for all j > 0, 

is an m-dependent sequence of vectors. We shall now show that {t;*/™\ j > 0} satisfy 
conditions similar to d?]) and ^ and then apply Step 2. Using (j98p and ([2]) in the equation 



'(A) = (2vr)-i/2^r;,,,-(t)J^(t/(m7,))e- 



we get that 

</"^(A) = ^J ^ K{mi)vl^{\ - 7-^e) . (99) 

It follows from Condition ([7]) that there exists a constant C > such that for all i > and 
AG [0,7r), 

|<,(A)| < C + |7,A - 7iA.,,|)-^ . (100) 
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Using ([To]) and the (27r)-periodicity of v* j, we can express (jlOOp using the symmetric (27r7j)- 
periodic function /i27r7j.7jAi ^ (0 defined in Lemma [2] and equal to (1 + — for 
< ^ < vr7j. With ^ = jjX, one gets 

|<,(A)| < C /i2.7„7.A,.,(7.-A) , j > 0, A G M . 

Let g{t) = m\K(mt)\ and observe that \\g\\i = \\K\\i < oo and 

g{t) < Co m < co\t\'^''^ for aU jtj > 1 and m > 1 , 

where cq is a positive constant such that < co|u|~'^~"'^ for \u\ > 1. Applying these 

bounds to (f99]) gives 



1/2 



5(0 ^27r7,,7,A.,,(7iA-0 



Applying Lemma [2] to this convolution, we get 

^•>\A)|<C7 7j/' (l+7,|A-A,,|)-^ 

for different constants C depending neither on m > 1, j > nor on A € [0, vr). One has 
therefore the following version of ([7]) for v* uniform in m > 1: 



sup sup sup 7"^/^|</™Ha)|(1 + 7j|A- A,j|)^ < oo 

j>0 m>l Ae[0,7r) 

To get a version of dHj) for v*j^'^\ observe that, by ([99]) . we have 



(101) 



m 
2^ 



de 



Condition ([7]) implies that the term in brackets is bounded independently of ^ and j and 
hence by ^ and dominated convergence, one has 



Iim77'/M/'"^(77'A + A,,) 



(m) 



(A) for ah A G 



where 



^(^o<oo(A-e)d^ 



(102) 



(103) 



Note that v* ^^^\X) is an approximating sequence of v*^{\) in the sense that, since v*^ is 
bounded (by (|14|) ) and continuous (by hypothesis), and since (27r)~^ -ft'(n) du = K{Q) = 
1, for all A e M, 



<oo^^'"^(A) 



27r 



K{u)vl^{X-u/m)du^vl^{\) as ?n ^ oo . (104) 
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Relations (jlOip and (jl02p are the corresponding versions of Conditions ([7]) and ([8|) for 
V* and since, we are in the m-dependent case, we may apply the result proved in 
Step 2, and obtain, as j oo, 

^nj — 1 J y{m) 2 

-1/2 



Z^fc=0 i-^lj.fc ~ J J" 



^nj-lr„(m)2 |- ^m;z1^ 

where F^"^) is the covariance matrix with entries 



,(m)2i 



(105) 



dA, l<i,i'<iV, 



(106) 



where w*^,^"^^ is the equivalent of w*^, in (fT6 



, AG 



To obtain the corresponding result for the {Z^ jk} sequence, we apply [3|, Theo- 
rem 3.2] as follows. We show that 

lim Vf!) =T,^i, l<i,i' <N (107) 



and, for alH = 1, . . . , A^, 



lim lim Var | n, 

m— *oo j— >oo 



-1/2 



E - E = . (108) 

fc=0 fc=o / 

Relation (jl07p says that the RHS of (jlOSh converges in distribution to the RHS of ([23]) 
and Relation (jlOSp says that the LHS of (llOSp is a good approximation to the LHS ()23|) by 
choosing m arbitrary large. 

To prove (jlOSp . it is sufficient to establish the following equalities for all i = 1, . . . , A'^. 



lim lim Gov I n ■ ^^"^ Z." 



Am)2 -1/2 



m— >oo »oo 



k ' J 



k=0 



k=0 



lim lim Var I n ■ ^''^ Z^^l'^ 1 = lim Var I n ■ ""^^^ ■ 



,{m)2 



1/2 



m — >oo J — ^oo 



k=0 



k=0 



Using Relation (jTSp of Corollary [H the limits as j — > oo (an hence nj oo) in the previous 
display are, respectively, 

2 



-p(m,oo) def „ „ 



* (m,oo) 



(A + 2p7r) 



dA, rj^) and Ti 
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where T^™^ is defined in ()106p and Fj^j in p4p and 

1 



2 L 



Hence to prove (|107p and (jlOSp . it is sufficient to show that 



lim rS™'°°)= Um T^=T,,,, 

m— >oo m— >oo 



= 1,...,N . 



Observe first that Relations (jlOip and (jl02p imply 



sup sup 

■m>l 



*J'")(A) (l + |A|)^<oo 



A G 



(109) 



(110) 



which is the uniform version of ((111). Eq. (1109P now follows bfrom (I104p . (jlOlD and (jllOp . 
and dominated convergence. □ 



7. Proof of Theorem [2] 

The following proposition is the key point for proving Theorem [2] since it shows how 
Condition ([7]) in Theorem [T] can be recovered for an approximation of the sample mean 
square, when using the alternative Condition ()26p . Condition (I27p in Theorem [2] can then 
be used to control the sharpness of the approximation. 

Proposition 4. Let {Zij^fc, i = 1,2, j > 0,k £ Z} be an array of 1 -dimensional decimated 
linear processes as defined by (Cp. Assume that {vij{t), j > 0, t € Z} satisfies (26\) for 5 > 
1/2, a sequence (Aij) taking its values in [0,7r) and some e > 0. Then there exists an array 
{vij{t), j > 0, t € Z} whose Fourier series coincide with those of {vij{t), j > 0, t € Z} 
in e -neighborhoods of the frequencies {Aij, j > 0} and satisfying that is, such that 



''^ijW — f'^'"' A € (— 7r,7r) such that |A — Aijl < e 

sup sup 7"^''Vt,-(A)|(l+7j|A-Aij|)^ < cx) , 

i>0AG[0,7r) 



(111) 

(112) 



and the following approximation holds. 



n. 



-1/2 



Y,{zi,,k-nzij,k\ 



k=0 



k=0 



where 



Zi,j,k = ^vijiljk - t)6 
and, for some positive constant C not depending on j. 



E[\R,\]<C 



1/2 1/2 

+ 1- 



:il3) 



:ii4) 



:ii5) 
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where 



Ij = / 1(|A - Ai,oo| > e) \vlj{X)\' dX . (116) 
Jo 

Proof. Let Lq = [-Ai,oo - -Ai,oo + e] U [Ai,oo - Ai,oo + e]- We write 

where vlj{X) = 1lo(A)w^j(A) so that pil|) holds. We define vij, vij accordingly, so that 
vij{t) = vij{t) + vij{t) and, since j and vl ^ are in L^(— 7r,7r), vij and i^ij are in 
Hence = Zi,j,k + -2^1,^,^ with Zij^k defined by (jll4p and 

Moreover, by (j26p and the definition of ujj'j. Condition (jll2p holds. 

We now show that the remainder Rj defined by (|113p satisfies (jllSp . Observe that Zij^f^ 
and Zi j f; are centered and, since 



E[Zij- fcZij-fc] = / t;* .(A)fi*^^.(A)dA = , 

J — TT 

uncorrelated. Thus we get E[Z^^. ^] = E[Z^^ J +E[Z^^. ^] and hence the remainder Rj 
defined by (fTTSD is 

Rj = Pj + 2Qj with (118) 



Pj=^j'^'T.i^hl^-^^^hk]} and = n//' ^ Zi,,, fcZi,,, . (119) 



k=0 k=0 

We have 

E[|P,|] < 27^fE[Zl^,] = 2nfY.^,/t) = Anf /, . (120) 

by the Parseval Theorem and the definitions of j and Ij . 

Using that Zij^k and -^ij,fc' are centered and uncorrelated, we have, using a standard 
formula for cumulants of products, for all k, k' € Z, 



Gov yZij^kZij^k, Zij^k'Zij^k'j = Gov yZij^k, Zij^k'j Gov yZij^k, Zij^k 

+cum (^Zi j k, Zij^k, ^i,j,k'i -^i,i,fc' 

Hence, Var (Qj) = Aj + Bj where 

rij — 1 iij — l 

= '^7^ X] X] (^i,i,fc>^i,j,fc') Gov (^Zij- fc, ^ij,fc' 

fc=0 A;'=0 
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k=0 k'=0 

Denote by fj and fj the respective spectral densities of the weakly stationary processes ^ij,. 
and Zij^,. Replacing the covariances in the defintion of Aj by their respective expressions as 
Fourier coefficients of the spectral density, e.g. Gov ^Zij^fc, Zij^fc/^ = e'^*'"'^ ^'^/j(A)dA, 
we get 



Ai = n • 



U\)I,(X') 



fc=0 



ifc(A+A') 



dAdV , 



which implies that 



0<^j<27r sup /j(A) X / /j(A')dA', 

AG(-7r,7r) 



where, in the last inequality, we used that, for any A', Y1^=q ^ 
Observe that, by definition of Zi j j. in (jll4p . 



i-1 ifc(A+A') 



:i2ii 



dA = 27rno 



Gov Zij-o,^ij,fc = (2^) / ^^i%(A) 



Using Lemma [T] with the (27r)-periodic function g{\) = v\ ^{X) 



i^ii^^^ we get 



Gov(Zi,,-o,Zi,,,fc) = (2vr)77^^" |^(77'(A + 2p7r)) 



e''^^ dA . 



Hence we have 



7,-1 

U\) = (27r)-S7i |^5(7-'(A + 2pvr)) 

Using ()112|) . since 17^""*^ (A + 2p7r)| < vr for p = 0, . . . , 7j — 1 and A e (— 7jvr, — 7jvr + 2ir), we 
get 



7,-1 

fj{X) <C ^ (1 + ||A + 2pTT\ - 7jAij|)-2^ A G (-7j^, -7^-^ + 27r) 

Using (j42p in Lemma [3] and that fj is (27r)-periodic, we obtain 

sup sup /j(A) < 00 . 

j>0 Ag(-7r,7r) 

Moreover, we have 

r fj{X')dX' = Var (Zi,,,o) = E^L'W = ^(^) ' ^A 



:i22) 



(123) 
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by the Parseval Theorem. Hence by (|12ip . there is a positive constant C such that 



\A,\ < C 



dA, 



(124) 



We now consider Bj. Using (A-2) and the definitions of ^ij.fc and in (|114p and (|117p . 
we have 



cum 



(^Zij^k, Zij^k, Zij^k' , = ^ vij i^jk- t)vij i^jk- t)vij {jjk' - t)vij {-fjk' - 1) 



Hence 



\Bj\ < ft:4 ^ \vij{t)vij{t)vij{t + -fjT)vij{t + -fjT) 

<Ki\y2 \^i,i{t)vi,j{'t)\ 



M / ki*,,(A) dAx / vl.{X) 



dA 



(125) 



By definition of v* ^ and (jl22p , we have 



vl .(A) dA = Var Zi,^- o = / /,(A) dA < C 



and hence 



B,<C 



dA 



(126) 



where C denotes a positive constant not depending on j. The bounds (I124p and (I126p and 
the definition of t;,* • yield 

E [Q2] = Var(Qj) = + By <2C j i^(A) ^ dA 

< 4C / 1 (|A - Ai,oo| > e) hi%(A)p dA = 4C 

by the definitions of ^ and Ij. This, with (jllSp . (|120p and Jensen's inequaUty yields (jllSp . 
which concludes the proof. □ 



Proof of Theorem\^ The general case can easily be adapted from the case A'^ = 1, which we 
assume here. We apply Proposition [H It follows from (jllip and (jll2p that the assumptions 
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of Theorem [1] are verified for Zij^k and we obtain 



n, 



-1/2 



k=0 
P 



It follows from i^, that Rj — >0 as j oo,. Hence (fml) yields the CLT ([23]), which 
concludes the proof. □ 

8. Proof of Theorem [3] 
We shall use the following lemmas. 



Lemma 7. Assume (H-1) and (H-2). For e > small enough, i/ 2/3 — 1 > 2, then 

2 



7i 



VF(7,A) /(A)dA-/(0) 



Oh] 



(127) 



Proof. Using (j3ip . for e > small enough, the left-hand side of (jl27p is at most 





Ty(7,A) 


' dA-/(0) 


+ c 






2 „ 

A^ dA 



(128) 



where C is a positive constant. To evaluate the first integral in (jl28p we write 
jT^ - fl^ - Using |(H-2)[ one gets 



7j 



Ty(7,A) 



dA = 1 + O ( 7; 



1-2/3 



The second integral in ()128p is bounded by 



7i 



VF(7,A) 



A^ dA < C'-f. 



where C = VF(A) A^ dA < 00 by[(H^ 



□ 



Lemma 8. Assume (H-2) and define Bj as in Then there is a positive constant C, 

such that, for all j > and A G (— tt, vr), 

2 



and, for any positive e, 



< C 



1. 



1-/3 



W{ijX) 



+ 7 



1-2/3 



(129) 
(130) 



sup l(|A|>e) \B,{X)\ = 0{i/^-^) 

Ae[0,7r) 



(131) 
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Proof. By (fT38]) . we have Bj{X) = 7j^^Ty(7jA) + Rj{X), where 



Using (H-2) , since /? > 1, we have 



sup |ii,(A)| < C 7j/' 5^(1 + {2p - IhjT,)-" = 0(t 



AG(-7r,7r) 



(132) 



p>0 



The bound (fT32]) gives (fT29D . which yield (fmH by using ] (H-2) [ The bound (fT30D follows 
from (fT29D and 



•Zll - ^21 I < 2|Z2| X kl - -^2! + \zi - Z2\ 
1/2; 



apphed with zi = Bj{X) and Z2 = 7^ W^(7jA) 



□ 



Lemma 9. Suppose that the assumptions of Theorem\^hold. Then for some arbitrary small 
e > 0, 



-1/2 

u 

j>0 |A|<£ 

-1/2,,* . -1 



sup sup7/'^|<,-(A)|(l + 7i|A|)'^< 



00 



with 



and, as j —> 00, 



Proof. We have, for all u € 



lim77^'X.(77 A) = t;]',^(A) /or a// A€ 



vl,^{X) = a*{mW. A€M, 
l(|A|>e) K,/A)pdA = 0(7j-'^) 

1 



(133) 
(134) 

(135) 

(136) 



1 

2^ 



^^^#(^^.(A + 2p^)) 

pel 



e™^ dA 



hence the term in brackets is the Fourier series of {W{-fJ^u), ueZ} and thus 

uez pGZ 
which is some times called the Poisson formula. Inserting this in (I139p . we obtain 

<,(A) = a*(A)i?,(A), AG(-7r,^), 

where Bj is a (27r)-periodic function defined by 

Bj{X) = ^jfw{j,{X + 2p7r)). 



(137) 
(138) 
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Applying (jl37p and (jl29p in Lemma [H (H-2) and that |a*(A)| is bounded in a neighbor- 
hood of the origin by (H-1) , we get, for some arbitrary small e > 0, if |A| < e, 



ij'^l.Wl < C\W{^,\)\ + O < C (1 + T.IAI)"'^ + O 

where C and the O-term do not depend on A, which implies (jl33|) . 
Applying (jl37p and (jl29p in Lemma [8l we have, as j — oo. 



^13 



□ 



where the limit holds by (H-1) This gives (I134p . 

Applying (fT37D . (fT3TD in Lemma[8]and /^^ |a*(A)|2dA < oo, we obtain (jl36p . 

Proof of TheoremlM By ([281) and ([30]), we have 
where 

= 77'^' E ^(77' - ^)) «W> ^ e z . 

i-GZ 

Thus {Zij^t, J ^ 0, t G Z} is an array of one-dimensional decimated linear processes as in 
Definition [TJ Moreover, the Fourier series ([2]) of vij{s) is 



vlAX) = r 



'1/2 * 



(A) E^(^7'(^))«"''"' A€(-vr,7r). 



(139) 



We let = 1 and Aij = Ai,oo = for all j > 0, which yields ([3]), (gD, (P and §i) 
m Condition [(C)| In view of Lemma [U Relation ([8]) in Condition |(C)| holds, as well as 
Relations (j26p and ()27p in Theorem [2] (recall that in that theorem. Relation (j26p replaces 
Relation d?]) in Condition |(C)P . Hence we may apply Theorem [2l 

We are now in a position to show first ()34p then (I36p . Applying ()19p . we have 

(140) 



^ [^h,o] = Var (Zi,j- 0) = / \vljiX)\' dX . 

J — n 

Using (jl37p . (|136p and then (jl30p . this gives for any e > small enough 
^ [Zho] = j\ W{X) B,{Xf dX + O (7)-'^' 



7j 



a*{X)W{^jX) 



dA + C 



w{^,x) dx + oU-^P) . 



In the last line, since |a*(A)|^ = /(A), by Lemma [3 the first term is /(O) -|- O and, 
by a change of variable, the second term is less than C77^||VF||i. Hence 
/3 > 2. 



follows since 
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The bound ()136p yields (f27I) under the condition 7- = o(n - ). Since nj is given 



by ()33p . the assumption ([35]) imphes that condition. Since (A-2) holds as well, we may 



apply Theorem [2] with j = n and obtain (j36p . □ 
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